Observation of polaritons in E^S^CaC^Og+z single crystals 



Sven-Olof Katterwe, Holger Motzkau, Andreas Rydh, and Vladimir M. Krasncrvj^] 

Department of Physics, Stockholm University, AlbaNova University Center, SE - 106 91 Stockholm, Sweden 

The Bi2Sr2CaCu2C>8+:r high-temperature superconductor represents a natural metamaterial, com- 
posed of metallic CuO bilayers sandwiched between ionic BiO planes. Each pair of CuO bilayers 
forms an atomic-scale Josephson junction. Here we employ the intrinsic Josephson effect for in 
situ generation and detection of electromagnetic waves in Bi2Sr2CaCu2 08+a; single crystals. We 
observe that electromagnetic waves form polaritons with several transverse optical phonons. This 
indicates the presence of unscreened polar response in cuprates, which may lead to strong electron- 
phonon interaction. Our technique can provide intense local sources of coherent, monochromatic 
phonon-polaritons with kW/cm 2 power densities. 



Propagation of light in a media is accompanied by 
electromagnetic (EM) interaction with charged particles. 
The interaction is strong in polar materials, such as ionic 
insulators and ferroelectrics. Resonances between light 
and EM-active collective modes in the media (e.g., op- 
tical phonons 1], or excitons [2l |3]) lead to formation 
of polaritons, half light - half matter particles, with a 
singular dielectric function 

e r H = eoo + E 2 ■ (1) 

It has poles at transverse wto and zeroes at longitudinal 
w lo optical frequencies. Here Sj and are oscillator 
strengths and damping parameters, respectively. Light 
does not propagate through the media at wto < oj < ^lo 
for the same reason that low frequency EM waves do not 
propagate through metals, namely due to the negative 
dielectric constant Re(e r ) < 0. For 7 — > the group 
velocity of light turns to zero at wto- The ability of po- 
laritons to slow down, or even stop and store photons has 
attracted significant attention in recent years because it 
may be used in optical metamaterials and facilitate 
manipulation of photons in optoelectronics [6] . 

Here we study EM-wave propagation in the layered 
high temperature superconductor Bi2Sr2CaCu208+z 
(Bi-2212). Superconductivity in cuprates occurs as a re- 
sult of doping of the parent ionic oxide insulator. Opti- 
cal spectroscopy has revealed an intermixture of metallic 
(Drude-type anomaly at uj = 0) and ionic [strong infrared 
(IR) phonons] properties in cuprates [?HS]. However, Bi- 
2212 is not a uniform doped Mott-insulator. It repre- 
sents a natural atomic superlattice, composed of metallic 
Cu02-Ca-Cu02 bilayers sandwiched between ionic SrO- 
2BiO-SrO planes. Most clearly this was demonstrated 
by observation of the intrinsic Josephson effect between 
CuO bilayers [TU1 - [T2"] . From an optical point of view, 
Bi-2212 represents a layered photonic crystal - 
a multilayer transmission line for EM waves. Therefore, 
low frequency transverse EM (TEM) waves can propa- 
gate along insulating BiO planes [TB] despite that Bi- 
2212 is metallic and superconducting. Importantly, the 
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intrinsic Josephson effect allows both in situ generation 
and detection of TEM waves in the crystal, i.e., the ma- 
terial as such can convert the dc-bias current into high 
frequency ac-electric field via the quantum mechanical 
ac- Josephson effect. 

We study small mesa structures, fabricated on top 
of Bi-2212 single crystals using micro/nano- fabrication 
techniques. We use two batches of crystals: pure, 
slightly underdoped Bi-2212 with critical tempera- 
ture T c ~ 82 K, and lead-doped, slightly overdoped 
Bi 2 _ a Pb y Sr 2 CaCu 2 8+2; with T c ~ 87-91 K. The most 
noticeable difference between them is in the c-axis crit- 
ical current density J C (4.2K) ~ 10 3 and 10 4 A/cm 2 for 
pure and lead-doped mesas, respectively. Measurements 
were performed in a gas-flow 4 He cryostat in a tempera- 
ture range down to 1.6 K and magnetic field H up to 17 T. 
Samples were mounted on a rotatable sample holder. Ac- 
curate alignment of the field parallel to CuO planes was 
needed to avoid intrusion of Abrikosov vortices in intense 
magnetic fields. For alignment we used the procedure de- 
veloped in Ref. [12] . which facilitates alignment with an 
accuracy better than 0.02°. Details of sample fabrication 
and measurements can be found in Refs. |TJ] and |16j . 

The inset of Fig. [I] shows current- voltage (I-V) char- 
acteristics of a lead-doped mesa (Ml, 1.3 x 1.9 /Ltm 2 ). Pe- 
riodic branches are due to one-by-one switching of junc- 
tions from the superconducting to the quasiparticle (QP) 
state [lOl [17] . Additional small sub-branches are seen at 
every QP branch. They represent resonances between the 
ac- Josephson oscillations in intrinsic junctions and c-axis 
LO phonons [TM2"T] . Figure [I] shows temperature depen- 
dence of the first QP branch for the same mesa. Unlike 
the superconducting gap and QP resistance, which de- 
crease with increasing T |17] , the voltages of phonon res- 
onances remain independent of T and correspond to fre- 
quencies of IR and Raman active c-axis phonons [7ll8ll22j. 
Their identification is provided in Table [I] 

The amplitude of phonon resonances is proportional to 
Jc [2] ■ Since J c for our lead-doped crystals is about ten 
times larger than for pure Bi-2212 crystals, phonon res- 
onances in the lead-doped mesas are greatly enhanced. 
For analysis of resonances in pure Bi-2212 we have to 
look on dl/dV characteristics [20]. Similarly, resonances 
are strongly suppressed by a moderate c-axis magnetic 
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FIG. 1. Temperature dependence of the first QP branch for 
the lead-doped mesa Ml in zero magnetic field. It is seen 
that phonon frequencies are T- independent. The inset shows 
QP branches in the I-V of the same mesa at T = 1.6 K. 
Resonances with LO phonons are seen at every QP branch. 



field. This is caused by penetration of Abrikosov pan- 
cake vortices in CuO planes, which suppress supercon- 
ductivity |23] and J c . In contrast, in-plane magnetic field 
up to 17 T is not affecting superconductivity and J c at 
T <C T c due to the extremely large anisotropy of Bi-2212. 
The apparent Fraunhofer modulation of the c-axis criti- 
cal current in this case 12J is caused by interference of 
opposite supercurrents, circulating in Josephson fluxons, 
with unsuppressed amplitude J c . 

At sufficiently large magnetic field, fluxons form a reg- 
ular lattice in the mesa [12) . The number of fluxons per 
junction is given by the ratio $/$o> where <£> = HLs is 
the flux per junction, L is the junction length, s ~ 1.5 nm 
is the stacking periodicity of intrinsic junctions, and 
$o is the flux quantum. The c-axis bias current ex- 
erts a Lorentz force and causes unidirectional fluxon mo- 
tion. This leads to appearance of an additional flux-flow 
branch in the I-V curve, with the voltage per junction 
Vff = uffHs, where mff is the fluxon velocity. Periodic 
collision of fluxons with the edge of the mesa leads to 
flux- flow oscillations [H] . Produced EM waves are partly 
emitted and partly reflected back into the junctions. The 
latter start to travel as TEM waves in the multilayered 
transmission line formed by the mesa, get reflected at the 
opposite edge, and continue bouncing until decay. 

Figure [2] shows flux-flow part of the I-V for a lead- 
doped mesa (M2, 1.2 x 2^m 2 ), measured during continu- 
ous sweep of the in-plane magnetic field from 4.2 to 4.4 T. 
The flux-flow voltage grows gradually with H . Bounc- 
ing TEM waves can be explicitly visualized at geometri- 
cal resonance conditions, corresponding to formation of 
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TABLE I. Identification of LO phonon resonances at zero field, 
marked in Fig.l (following Ref. 151). 
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FIG. 2. Flux-flow part of the I-V for the lead-doped mesa M2, 
measured during continuous sweep of the in-plane magnetic 
field from 4.2 to 4.4 T. The inset shows a series of Fiske 
steps, caused by formation of standing waves inside the mesa. 
It demonstrates the Josephson flux-flow mechanism of TEM 
wave generation in our experiment. 



standing waves in the mesa. They are seen as a series of 
Fiske steps on the flux-flow branch [16]. We note that 
those TEM waves are photon (or surface plasmon [5]) 
like, because they have a linear dispersion relation at low 
frequencies. This leads to equidistance of Fiske steps |16| . 
The inset in Fig. [^demonstrates the Fiske step sequence, 
corresponding to out-of-phase oscillations in neighboring 
junctions with one node per junction |16j . The corre- 
sponding slowest speed of light cn ~ 3.9 x 10 5 m/s is 
almost one thousand times slower than the speed of light 
in vacuum - not because of the extraordinary large e r , 
but due to an extraordinary large kinetic inductance of 
the atomically thin intrinsic junctions. 

The strongest geometrical resonance occurs at the 
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FIG. 3. A set of I-V curves for the lead-doped mesa M3 
at different in-plane magnetic fields and T = L6K. The 
curves are offset vertically for clarity. Development of the 
velocity-matching flux-flow resonance (circles) and higher 
flux-flow/phonon resonances (squares, rhombuses and trian- 
gles) is seen. With increasing H voltages of flux-flow-phonon 
resonances first increase and then saturate when the ac- 
Josephson frequency is approaching the next TO phonon fre- 
quency (marked by vertical lines). 



velocity-matching condition, when uff is equal to the 
phase velocity of TEM waves [TB] . In this case waves are 
synchronized with the fluxon chain and the wave length 
is equal to the separation between fluxons. The corre- 
sponding TEM wave number is 



k = 2irHs/^ a - 



(2) 



It can be easily tuned by changing H. In experiment 
the velocity matching resonance is seen as a sharp up- 
turn (maximum in dl/dV) at the top of the flux-flow 
branch. The corresponding frequency is given by the ac- 
Josephson relation 



2ttF; f /$o, 



(3) 



where Vp* F is the measured maximum flux-flow voltage 
per junction. Eqs. (2[3) provide the basis for analysis of 
the dispersion relation in this work. 

Figure [3] shows I vs. V per junction at different 
in-plane magnetic fields for a lead-doped mesa (M3, 
0.9 x 1.3 fJ.ro 2 ) with N = 11 junctions. The velocity- 
matching flux-flow resonance is marked by circles. It is 



moving to higher voltages with increasing H and even- 
tually saturates upon approaching the voltage of phonon 
#2 in Table |TJ Above it we observe a series of additional, 
remarkably strong resonances. The first three are marked 
by squares, rhombuses and triangles, respectively. They 
develop gradually out of one of the LO phonon resonances 
at H = 0, grow in amplitude and move to higher volt- 
ages with increasing field. Finally, they saturate at high 
fields when the voltage is approaching the next phonon 
resonance, as indicated by vertical lines. Those reso- 
nances can not be associated with bare LO phonons, 
which arc field independent, but represent a combina- 
tion of flux-flow and phonon resonances. Each series have 
N sub-branches with voltages V 3 ■ — (N — j)VpF + jVn, 
(j = 1,2,... N) due to one-by-one switching of N junc- 
tions from the flux- flow Vff to the resonance state Vr- 
N = 11 sub-branches in each series can be distinguished 
at H = 17 T in Fig. [3} The last sub-branch in each series 
represents the flux-flow/phonon resonance frequency. 

Figure [4] presents our main result: voltage per junc- 
tion of flux-flow/phonon resonances vs. H, which are 
translated into the dispersion relation uj vs. k via 
Eqs. (2p|. Figure Qa) shows the dispersion relation for 
the mesa M3. The lowest dispersion branch corresponds 
to the regular flux-flow branch, indicated by circles in 
Fig. [3j It is roughly linear at low k, corresponding to a 
constant speed of light ~ 4 x 10 5 m/s close to the out-of- 
phase TEM velocity, obtained independently from analy- 
sis of Fiske steps [16]. At larger k (H > 12 T) we observe 
a clear saturation with the group velocity duj/dk — > 
when the flux-flow frequency is approaching the lowest 
Raman active phonon #2 in Table [I] Higher dispersion 
branches correspond to flux-flow/phonon resonances, in- 
dicated in Fig. [3] They start flatly from a certain LO 
phonon frequency at k = 0, grow at intermediate k and 
saturate again at large fc, when the frequency approaches 
the next phonon (indicated by horizontal lines). Fig- 
ure (4^b) demonstrates a similar behavior of mesa Ml, 
which is slightly more overdoped and has about 30% 
larger J c . Figure |4|c) shows the dispersion relation ob- 
tained for a pure Bi-2212 mesa. Since phonon resonances 
are much weaker here due to smaller J c , we are able to 
trace interaction only with the two lowest phonons. 

The obtained results are in excellent agreement with 
theoretical prediction by Preis etal. [53] on the mutual 
influence of flux-flow and phonon resonances in Joseph- 
son junctions with ionic barriers. Interaction of oscil- 
lating electric field with IR and Raman active phonons 
leads to formation of polaritons, which slow down the 
group velocity of light to zero duj/dk — > at uj — > ujto 
and prohibit propagation of TEM waves and flux-flow at 
voltages <&oWTo/27r < Vff < 3>oWlo/27t. This leads to 
the observed stratification of I-V characteristics in Fig. [3] 
Our data is well described by the polaritonic dispersion 
relation of Eq. 0, as shown in Fig.Qd). 

Observation of strong polariton-like interaction of 
TEM waves with IR and Raman active c-axis phonons 
highlight the unusual natural metamaterial-type elec- 
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FIG. 4. Obtained TEM-wave dispersion relations for (a) 
mesa M3, (b) mesa Ml, and (c) a pure Bi-2212 mesa. It 
is seen that the dispersion relation is linear at low frequencies 
(constant speed of light), but saturate (zero group velocity) 
whenever the frequency is approaching the phonon frequency, 
indicating formation of polaritons. Panel (d) shows simulated 
polaritonic dispersion relation from Eq. |l]). Horizontal lines 
represent all observed LO frequencies. 



tronic structure of cuprates, in which two-dimensional 
metallic CuO layers are sandwiched between ionic layers. 
The combination of long range Coulomb interaction in 
two-dimensional metallic layers with polar media around 
them, may lead to an unusually strong electron-phonon 
interaction |25H29j . This is import for understanding of 
the puzzling phenomenon of high temperature supercon- 
ductivity. 



The maximum power, pumped into phonon-polariton 
emission can be explicitly estimated from the amplitude 
of flux-flow/phonon resonances in Fig. 3: P = A/Vr ~ 
10 fiW, where AI ~ 100 //A is the maximum resonance 
amplitude and Vr ~ 100 mV is the total voltage at the 
resonance. This corresponds to kW/cm 2 emission power 
density. It can be further enhanced in a coherent (cas- 
cade) manner using mesas with larger amount of stacked 
junctions [30] . Thus, our results demonstrate that volt- 
age biased Bi-2212 mesas act as compact and intense 
source of coherent, monochromatic phonon-polaritons (or 
surface plasmons) with power densities ~kW/cm 2 . Such 
local surface-plasmon sources are on demand in a rapidly 
growing area of optoelectronics |3HS]- The discovered 
phenomenon may also be used for creation of a THz 
phonon laser [30] . 
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